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Abstract 

We evaluate the fermionic Casimir effect associated with a massive fermion confined 
within a planar d + 1 dimensional slab-bag, on which MIT bag model boundary 
conditions of standard type, along a single spatial direction, are imposed. A simple 
and effective method for adding up the zero-point energy eigenvalues, corresponding 
to a quantum field under the influence of arbitrary boundary conditions, imposed 
on the field on flat surfaces perpendicular to a chosen spatial direction, is proposed. 
Using this procedure, an analytic result is obtained, from which small and large 
fermion mass limits, valid for an arbitrary number of dimensions, are derived. They 
match some known results in particular cases. The method can be easily extended 
to other configurations. 
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1 Introduction 



The macroscopically observable vacuum energy shift, associated with a quantum field, 
is the regularised difference between the vacuum expectation value of the correspond- 
ing Hamiltonian with and without the external conditions demanded by the particular 
physical situation at hand. At the one-loop level, when the external conditions are rep- 
resented by boundary conditions, this leads to the usual Casimir effect [l|] - see Ref. @ 
for an updated review on the theoretical and experimental aspects of this remarkable 
phenomenon. 

The physical fact that only bound states of quarks are experimentally observed, led 
to the model-dependent idea of the total confinement of the quark fields, through the 
so-called bag model. In general, the bag is modeled by the interior of a spherical shell 
within which the quarks and gluon fields are confined. There are several important cal- 
culations concerning the vacuum energy associated with the spherical bag model [§], as 
well as results concerning the influence of different boundary conditions imposed on the 
confined fermionic and gluonic fields ||]. Planar bag models with standard MIT boundary 
conditions (no fermionic currents through the surface of the bag) have been also consid- 
ered in the literature, in particular, the case of a massless fermionic field was considered 
in Ref. || and its extension to d + 1 dimensional was discussed in Ref. ||. The massive 
case in three spatial dimensions was considered in Ref. 0. 

In the context of the evaluation of the zero point energies associated with confined 
quantum fields, some configurations, which depend on the nature of the quantum field, 
the type of spacetime manifold and its dimensionality, and the specific boundary condition 
imposed on the quantum field on certain surfaces, lead to relatively simple spectra, easy 
to deal with, but others lead to rather complex ones. The heart of the matter in these 
calculations is the (physically meaningful) evaluation of the spectral sum that results 
at the one-loop level from the definition of the Casimir energy. This evaluation requires 
regularization and renormalization, and recipes for accomplishing this task range from the 
relatively simple cutoff method, employed by Casimir himself |IJ, to a number of powerful 
and elegant generalised zeta function techniques ||. Contour integral representations of 
spectral sums are a great improvement in the techniques of evaluating zero-point energies; 
they are especially useful when the spectra are not simple, and have been employed 
profusely before ||. 

Here, our aim is twofold. First, with the same objectives as in Refs. [|] in mind, we 
want to generalise to n + 1 dimensions a simple and effective way of evaluating the one- 
loop vacuum energy under external conditions, based on well-known theorems of complex 
analysis, namely, the Cauchy integral formula and the Mittag-Leffler expansion theorem 
in one of its simplest versions [11]. The method we will present, although very simple, 
is of sufficient generality so as to be successfully applied to a variety of cases. Second, we 
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want to apply our basic formula to the case of the n + 1 dimensional slab-bag, in which 
the quark field is constrained into a region bounded by two infinite hyperslabs separated 
by a fixed distance £. MIT boundary conditions (see Ref. and references therein) 
will be imposed on the quantum field. There are no fermionic field oscillations outside 
the interior region of the slab-bag, a situation that must be compared with the one in 
which the boundary conditions split the spacetime in several regions, and the zero-point 
oscillations of the relevant quantum fields are present in all of them. This example is one 
of the simplest one can think of that leads already to a relatively complex spectrum, and 
serves as a convenient test of the calculational tool that we present here. The paper is 
divided as follows. In Sect. 2 we introduce the method, and get the basic formula. In 
Sect. 3, we apply our main result to the problem in hand, and obtain its exact solution 
in terms of Whittaker functions. Finally, we consider some particular limits and check 
coincidence with some other results. The last section is devoted to final remarks. We 
employ natural units {% — c — 1). 

2 The unregularised Casimir energy and a very simple 
sum formula 

Consider a quantum field living in a d + 1 dimensional ultrastatic flat spacetime under 
boundary conditions imposed on two hyperplanes of area L d ~ l kept a distance t apart. 
Hence the motion along one of the spatial directions, say, the OAVaxis, is restricted. 
Suppose that the condition L 3> £ holds. At the one loop-level, the (unregularised) 
Casimir energy is given by 

^0 (d) =<x{d)—jY l ^£^n, (1) 

where a (d) is a dimensionless factor that counts the number of internal degrees of freedom 
of the quantum field under consideration, p± = J pi + p\ + ' ' ' + Pd-x-, an d 




where A is the n-th real root of the transcendental equation determined by the boundary 
conditions, £ is a chararacteristic length along the OX& direction, and m is the mass of an 
excitation of the quantum field. A simple integral representation of J2 n can be written 
if we make use of Cauchy's integral formula. In fact, it is easily seen that 
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where in principle T is restricted to be a Jordan curve on the g-complex plane with ^sq > 0, 
consisting of a semicircle of infinitely large radius, whose diameter is the entire real axis. 
Taking (|3|) into ([!]) we obtain 

E ° {d) = ~ a {d) — J (2^1 fr^? (4) 

In order to proceed we must be able to perform (in a physically meaningful way, as 
mentioned above) the discrete sum in (||). 

Consider a complex function G(z) of a single complex variable z, symmetrical on the 
real axis and such that its roots are simple, non-zero and symmetrical with respect to the 
origin of the complex plane. The assumption that the origin is not a root of G(z) is not a 
restrictive one, because if z = happens to be a root of G(z), we can always divide G(z) 
by some convenient power of z in order to eliminate zero from the set of roots, without 
introducing any new singularity. Let us order and count the roots of G(z) in such a way 
that 

A n = -A_„, n = ±l,±2,.... (5) 
Now define the following meromorphic function 

oo / -i 

JW = E — P (6) 

n=-oo z lA n 

where the prime indicates that the term corresponding to n = must be omitted from 
the sum. The following properties of J(z) are straightforward: (i) J(z) has first order 
poles which are determined by the roots of G{iz): (ii) the corresponding residua are all 
equal to one. Taking into account property we see that J(z) can be rewritten as 

1 / ~ ' 1 ~ ' 1 \ 




z + i\ n 



Let us consider now the function K(z) := G{iz). We can state that 



(7) 



J{z) - (8) 

where the prime stands for 'derivative with respect to z\ In fact, the function K'(z)/K(z) 
has the same simple poles as the originally defined J(z), and also the same residuum at 
each pole. Hence, we can invoke the Mittag-Leffler theorem, to conclude that Eq. <M) is 
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true up to an entire function which does not contribute to the evaluation of the vacuum 
energy. It then follows that we can write 

K'(z) d , rr „ ~ 2z 

In order to make use of Eq. @, we first link the dimensionless complex variable z to 
the auxiliary complex momentum variable q, through the relation 

f + & n = (10) 

hence 



z = z (q,p±) = t\jq 2 +p 2 ± + m 2 , (11) 

and 

00 2a 2 C 2 o 2 00 2z 

V— - = — — V (12) 

^ a 2 + Q 2 z ^ z 2 + \ 2 ' y ' 

n =l <i ^"n z n=l Z ^ A n 

Notice that A = is explicitly omitted on the lhs of Eq. (|T^). Changing variables (d/dz = 
(z/£ 2 q)d/dq), we obtain for the unregularised Casimir energy the following expression 



ft M = ^ M ^ / g 8 | (13) 
which can be integrated by parts to yield 



e oW = - a(d) ^/g£i/ r gl {? i„ g[ ^) 1} 

Notice that the integration is now performed on an open curve which lies on the Riemann 
surface of the integrand the projection of which on the g-complex plane is the curve T. 
This curve begins at a point in the second quadrant very close to the imaginary axis and 
very far from the origin, describes an arc of a circle counterclockwise until it meets the 
real axis, then proceeds to the right along the real axis towards a point very far from 
the origin, where it starts to describe another arc of circle counterclockwise until it meets 
again the imaginary axis (Fig. 1). 

The first term on the rhs of (|i~4"D contributes with a phase which cancels out with a 
phase coming from the second term, being the final result real. The second term still 
needs regularization and in order to accomplish this we split the function K (z) in two 
parts 

K(z) = K 1 (z) + K 2 (z), (15) 
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Figure 1: The curve T begins at a point in the second quadrant very close to the imaginary 
axis and very far from the origin, describes an arc of a circle counterclockwise, then proceeds to 
the right along the real axis towards a point very far from the origin, where it starts to describe 
another arc until it (almost) meets again the imaginary axis. The starting and end points are 
symmetrical with respect to the imaginary axis. T\ is the part of T that goes from the starting 
point to the origin and is the part from the origin to the finishing point. 



with the following properties: (i) all terms whose integrals diverge when Re z > are in 
K 1 and (ii) all terms whose integrals diverge when Rez < are in K 2 , (iii) the symmetry 



K x (z) = K 2 {-z) , 



(16) 



holds. Then, we can write 



/ := 



dq 
r 2tt 



log [K(z)\ 



dq 



1 + 



Ko(z) 



dq 



1 + 



K x (z 
Ki(z) 



(17) 



where the integrating path has been divided into two parts: Ti is the part of T that goes 
from the starting point to the origin and T 2 is the other part, from the origin to the end 
point (Fig. 1). The integrals over logi^ (z) , j = 1,2, are evaluated along the different 
parts, Tj, of the curve T and, owing to the symmetry of the construction, their sum is 
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equal to the integral of \ogKi along the total path T; this yields a phase that, as remarked 
above, cancels the phase coming from the first term in Eq. ([HI) . The other two integrals 
are reduced to their contribution along the real axis, and we can write 



d-l 



E Q (d) = a (d) 



d d p 
(2vr) a 



log 



1 + 



Ki{z) 
Ko(z) 



(18) 



where all along the real axis z does not change its sign and is a function of p '■= {q '■= Pd, P±) 
and m given by z = +£\/p 2 + m 2 . Notice that Eq. flilf) gives the regularised confined 
vacuum energy. 



We can proceed a little bit further, still without specialising Eq. (|i~8|). First we make 
use of 



d d pf(p) 



d 

2tt2 



r (i) 



dp P d - i f ( P ) 



(19) 



and rescale the integration variable according to p — > xjl to obtain 



E H,d) = a (d) 



-d-l 



(2vr) d r 



I) £ d Jo 



dx x d 1 log 



1 + 



Ko(z) 



(20) 



where now z = +y/x 2 + /i 2 with /i := m£. Let us change the integration variable, accord 
ing toz^w 
by 



x 2 + /i 2 ) 1//2 . Then it follows that the regularised vacuum energy is given 



Eq (£, /i, d) — a (d) 



L 



d-l 



(2vr) d r 



dio (to 2 



2\ 2 



lot 



1 + 



(21) 



Equation ([21]) - or if one wishes Eq. (0) - is our main result. It is not hard to realize that 
this result holds for arbitrary boundary conditions imposed on the field on flat surfaces 
perpendicular to the OX^ axis. However, due to the subtraction of the zero root from 
the set of all roots, an extra multiplicative factor of 2 will be needed when dealing with 
topological periodic or antiperiodic conditions. The function K can be constructed from 
the boundary conditions, as we will see explicitly in the next section. 

3 Massive fermion field under MIT boundary condi- 
tions 



As an example of the usefulness of (|18D , let us apply it to a massive fermion field confined 
within an hypothetical hyperbag with MIT boundary conditions imposed on the field, on 
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the OX d axis. It can be shown that, for this relatively simple bag model, the eigenvalues 
of pd are determined by the roots of the function F(xd '■— Pdty, defined by 



F(pd) = ^sm(p d £) + p d £cos(p d ( 
Therefore we can choose G(uj) as 

sin (uj) 



G(uj) = ji- 



+ cos (uj) 



uj 



(22) 



(23) 



where we have divided F(uj) by uj because uj = is a root of F (uj). We can easily prove 
that for the case at hand the roots of G{u) are all real. Now, we construct K {uj) in 
accordance with 

K(uj) = G(iu) = /i Smh ( ^ + cosh (uj) . (24) 



UJ 



It follows that we can write 



and 



A'i M 



-11-^1. 

UJ 



K 2 (uj) = \(l + ^ 
2 V uj 



(25) 



(26) 



Notice that the symmetry given by Eq. (jig ) is obeyed as it should. For a fermionic 
quantum field in d + 1 dimensions, we have a (d) = —2xC (d), where the factor 2 
takes into account the existence of particle and antiparticle states, while C (d) counts the 
number of different spin states and is given by 2 ( - d ~ 1 ^ 2 1 for d odd, and by 2( d ~ 2 ^ 2 , for d 
even. Therefore the regularised vacuum energy of the fermion field is 



-d-i 



E (£, fi, d) = -C (d) 



2 rf-l 7r d/2 r f | 



duo ( 



2 2\ 2 

UJ — /! 1 



log 



UJ — u 9 , , 

1 H e~ 2ul 

uj + /i 



(27) 



Setting d = 3 we readily obtain the result given in Ref. (?]]. Equation (|27l) can be integrate 
by the following procedure. First we make use of 



io g (i +*) = -£; 



;-x) A 



(28) 



k=l 



which holds for — 1 < X < 1, and write Eq. ( PH| ) under the form 



E (£,fx,d) 



C(d) L 



d-l 



E 

fc=i 



-h (a», c^) 



(29) 
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where 



h (/•*, d) :-- 



duo uo (u + n) 2 1 k (uo — jj) 



In order to evaluate this integral, consider the auxiliary one 



I k (//, d, A) := 



duo (u + /i) 2 1 k (uo — fi) 2 



so that 



1 d 



-2ku}\ 



4 (/i, d) = -i^-^h 0", d,X = l). 



(30) 



(31) 



(32) 



To evaluate the auxiliary integral we make use of (c.f. formula 3.384.3 in [|i~2|l ): 

"(/3-//i) 



, / . n\2v-l / \2p-l -u 2 x (/^l + ^) +P 

da; (£ + /}) (x-/ii)^ e ^ 2 = ^ exp 



mi 



x r (2p) ^_ Pi ^ +p _ i (^i/i 2 + Pfr) , 



(33) 



where W Utfl (z) is the Whittaker function ||13|| . This result holds for /xi > 0, |Arg (/3 + /ii)| < 
Ti, Re /i 2 > and Re p > 0. The auxiliary integral then reads 



— — 1 / \ 

I k (/i, d, A) := -^^r ( ~ + *W_ fcl fci (4^A) . 
(2fcA) 2 V 2 / 



(34) 



From this result we easily obtain 



h d) 



r(| + fc)^ j 

4fc^ rfA 



A -2 (4//JfeA) 



A=l 



(35) 



Therefore the Casimir energy is, at last, 

c (d) L<*-y 



2 d+l 7l d/2 T ^ 

d 

X dA 



E 



\-i W_ k 4=1 (4/xJfeA) 



A=l 



(36) 



Analytically, this is as far as we can go and can indeed be considered a nice result. In 
order to get, however, a taste of the order of magnitude of the different terms involved, 
or to get a graphical plot of the Casimir energy, we must evaluate Eq. (|36|) numerically 
(which can be done very quickly) or consider some special limits. 
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3.1 The limit 



The limit \x <C 1 is determined by the behavior of the Whittaker function on the positive 
real axis for small values of the argument. For completeness, we state some facts concern- 
ing this function and relevant to our purposes in the Appendix. It follows then that, for 
/i < 1, the Casimir energy is given by 



E (i,fJ.,d) 
L d-i 



C{d) 


[2 d - l) 


T(d)C R {d+l) 


2 3d-l 7r d/2 r 


3) 


£ d 



(37) 



(38) 



Setting d = 3 we obtain the well known result due to Johnson ||, namely 

E (£, /i, d) _ 7tt 2 
TF 1 ~ ~2880£ 3 ' 

Equation ([37]) is also in agreement with the analysis of De Paola, Rodrigues and Svaiter 
for the d+1 dimensional slab-bag for massless fermions 0. This can be easily verified by 
setting d = D — 1 and making use of the duplication formula for the gamma function ||i~2| 



2? x ^ / 1 

T(2x) = ^^T(x)T[x + -), (39) 



with x — (D — 1) /2. For d = 3, in order to compare our result with one obtained by 
in this limit we take into account the second term of the small z approximation of the 
Whittaker function (see the Appendix). The result is then 

2880^ (1 " Bfl) ' (40) 

where B = 540 ( (3) /7tt 4 w 0.95 which is approximately half of the value obtained in [0]. 



3.2 The limit /i > 1 



Making use of the large |^| representation of the Whittaker function (see the Appendix), 
we obtain 



E {£,fi,d) 



-C(d) 



+ 1 



d 



Setting d = 3 and keeping only the leading term, we get 



- + 1 U" 2 + 2/i^ 1 



-2/i 



Eb(^,/i,d=3) 



3/i 1 / 2 



-2/x 



(41) 



(42) 



L 2 2 5 7T 3 / 2 £ 3 

in agreement with as soon as we take into account the particle and antiparticle states 
and thus multiply these authors' result by a factor of 2. 
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4 Conclusions 



In this paper we have derived a general regularised expression, Eq. (0), for the evaluation 
of the Casimir energy of a quantum field in a flat manifold under the influence of boundary 
conditions, imposed on the field on flat surfaces, or topological conditions constraining 
the motion along a particular spatial direction. Then we have applied the main result to 
the case of a massive fermionic quantum field confined by a planar n + 1 dimensional slab- 
bag with MIT boundary conditions. An exact result for the confined fermionic vacuum 
energy was obtained in terms of the derivative of the Whittaker function with respect to a 
dimensionless parameter. Approximate expressions for the vacuum energy were obtained 
in the limits of small and large values of the parameter // = £m. 



The authors have verified explicitly that the simple Eq. (|18|) works perfectly well in 
several other instances. For example, it also holds when the OX d direction is compactified 
by the imposition of topological conditions, periodic or antiperiodic, on the quantum field 
along that direction. In fact, the method can be extended and applied to cylindrical and 
spherical geometries embedded in d + 1 dimensional spacetimes and, possibly, to other 
more complex combinations of fields and geometries. 
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Appendix: The Whittaker function 

The Whittaker function is defined by pi 



W ua (z) := e~ z ' 2 z a+l ' 2 U Q + a - u, 1 + 2a; 



where U (a, b; z) is the confluent hypergeometric function of the second kind, 

r (a) L ' 



U (a, b; z) = f e~ zt e- 1 (1 + tf'^ 
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For the situations we are interested in, the small \z\ behavior of this function is given by 

U(a,b : z) w T ^~ 1 ' z 1 - b , Reb>2,b^2. 
r (a) 

As a consequence, for small \z\, the Whittaker function behaves as 

r ( 2<T ) f .-a+l/2 _ 1 -0-+3/2 



W va (z) » . v y N - ~z 



2 ' 



2 



For large values of \z\, the behavior of the Whittaker function is given by a series repre- 



sentation [12], fL3|, whose leading term is 

W ua (z) « e- 2/2 ^. 
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